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Abstract 

First, we investigate the static non-Abelian Kubo equation. We prove that 



CN . it does not possess finite energy solutions; thereby we establish that gauge 

^^ , theories do not support hard thermal solitons. This general result is verified 

by a numerical solution of the equations. A similar argument shows that 

o, 

'nJ" , "static" instantons are absent. In addition, we note that the static equations 

ON 

reproduce the expected screening of the non-Abelian electric field by a gauge 



Oh! invariant Debye mass m = gT J „^' . Second, we derive the non-Abelian 



Kubo equation from the composite effective action. This is achieved by show- 



k> i ing that the requirement of stationarity of the composite effective action is 

H : 

5t , equivalent, within a kinematical approximation scheme, to the condition of 

gauge invariance for the generating functional of hard thermal loops. 
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I. INTRODUCTION 

When it was realized |l[] that the gauge invariance condition on the generating func- 
tional T{A) for hard thermal loops in a gauge theory Q] (with or without ferniions) coincides 
with a similar requirement on the wave functional of Chern-Simons theory, one could use 
the known solution for the latter, non-thermal problem ||^ to give a construction of T{A) 
relevant in the former, thermal context. The expression for T{A) is non-local and not very 
explicit: T{A) can be presented either as a power series in the gauge field A [|I| [the OlA"^) 
contribution determines the hard thermal gauge field (and fermion) loop with n external 
gauge field lines] or as an explicit functional of path ordered variables Pexp / dx^ A^ [§]. 

More accessible is the expression for the induced current — ^j^ = — T" ^^^ , which enters 

H- (J. 

(high-temperature) response theory, in a non-Abelian generalization of Kubo's formula (in 
Minkowski space-time) 0: 

T" is an anti-hermitian representation of the Lie algebra, the gauge covariant derivative is 
defined as Dy = d^^ + g[A^, ], and m is the Debye mass determined by the matter content: 
in an SU{N) gauge theory at temperature T, with fermions in the representation T°, and 
Tr (T"'T^) = — ^5'^^ where Np counts the number of flavors, the Debye mass satisfies 

Henceforth, through Section |I|, we scale the gauge coupling constant to unity. The functional 
form of j^ can be given as [^ 

fix)=J^iQ't{a^{x)-A^ix)^+Q^_(a^ix)-A^ix)^Y (1.3) 

Here Q± are the light-like 4-vectors -y^i^, i^); with g^ = 1, A± are the light-like projections 
^± = Q± ^ni while a± are given by [^^ 

a+ = g-^d+g, a. = h'^ d^ h {d±=Q^df,) (1.4) 



when A± are parametrized as 

A+ = h-^d+h, A_=g-^d^g . (1.5) 

In other words, a± satisfy the equations 

a+a_ - d-A+ + [A+, a-] = , (1.6a) 

d+A_ - d.a+ + [a+, A_] = , (1.6b) 

whose solution can be presented as in (|1.4|) when A± are parameterized as in (|1.5|) — evidently 
g and h involve path ordered exponential integrals of A±. (Alternatively a± may be given 
by a power series in A^ |jl|.) Finally ( |1.3|) requires averaging over the directions of q. 

It is easy to verify that (|1.6| ) ensure covariant conservation of j'^. Moreover, gauge 
invariance is maintained: for ( p..l|) to be gauge covariant, it is necessary that j'^ transform 
gauge covariantly. That the expression in (|1.3| ) possesses this property is seen as follows. 
When A± transform by U^^ A± U + U^^d±U, Eqs. ( |1.4| ) - ( |1.6| ) show that a± transform 
similarly, hence the differences a±—A± transform covariantly. The manifest gauge covariance 
of ( |1 . 1[ ) ensures that m is a gauge invariant parameter; that it also has the interpretation of 
an electric (Debye) mass will be evident when we consider the static limit. 

It is of obvious interest to discuss solutions of (|1 . 1|) . In the Abelian, electrodynamical 
case this is easy to do, since ( pTBI) can be readily solved for a±, and the solutions of the 
linear problem are the well-known plasma waves . The non-linear problem of finding non- 
Abelian plasma waves is much more formidable. Also, one inquires whether the non-linear 
equations support sohton solutions, and (after an appropriate continuation to imaginary 
time) instanton solutions. [The time-dependent equation (|1 . 1|) in Minkowski space-time 
must be supplemented with boundary conditions, which are determined by the physical 
context. For example, response theory requires retarded boundary conditions, which in fact 
preclude deriving ( |1.1| ) variationally 0. Here we shall not be concerned with this issue.] 

Our paper concerns the following two topics. In Section |I|, we analyze ( |1.1D for static 
fields. It turns out that in the time-independent case (|1.6| ) can be solved for a± and (|1.1| ) is 



presented in closed form. We prove that the resulting equation does not possess finite-energy 
solutions, thereby establishing that gauge theories do not support hard thermal solitons. 
Also some negative conclusions about instantons are given. In Section |I1| we present an 
alternative derivation of ( |1.1| ), which relies on the composite effective action |]^, and makes 
use of approximations recently developed in an analysis of hard thermal loops based on the 
Schwinger-Dyson equations [§]. The Appendix presents a numerical analysis of the solutions 
to equation ( |1.1|) for SU{2), which supports the general result in Section |T|. 

II. STATIC RESPONSE 



When A± are time- independent, we seek solutions of ( |1.6| ) that are also time-independent. 
Acting on static fields, the derivatives d± become ±-4jg- V = ±9,-, and ( |1.6| ) may be written 
as the equations 

drA±±[A±,A±]=Q (2.1) 

for the unknowns A± = A± + a^. These are solved trivially by A± = 0, that is 

a^ = -v4± . (2.2) 

This solution is also the one that is deduced from the perturbative series expression for a±, 
when restricted to static A±. 

[A non trivial solution can be constructed with the help of representations similar to 
[Si). Upon defining in the static case 



y4+ = /i/ dr ho , 
A- = -Qo^drgo 

{ho and go involve path-ordered exponentials along the path r + gr), we find 

A+ = h^^ /+ ho , 

A- = go^ I- go , 



where I± are arbitrary Lie algebra elements, independent of r: q ■ V/± = 0. Since these 
solutions involve the arbitrary quantities J-t, and since they do not arise in the perturbative 
series, we do not consider them further and remain with the trivial solution ( |2.2| ), which 
corresponds to /± = 0.] 

From ( p^.2D it follows that the current for static fields is 

= -/^(Q'; + e)(Q^ + g^)^.(r). (2.3) 

With Q+ + Q_ = and (5+ + <5° = v^, we compute j^ = —2 5^^Ao. The response equations 
( |1.1D then become, in the static limit: 

DiE' + m^Ao = , (2.4a) 

e'^'^DjB'' = [Ao,E'] , (2.4b) 

where E^ = F'^^ and F^^ = —e^^^B^. Eqs. (|2.4| ) give clear indication that m plays the role 
of a gauge invariant, electric mass. The fact that the static current is linear in the vector 
potential implies the vanishing of hard thermal loops with more that two external gauge-field 
lines, and zero energy — a fact which can be checked from the relevant graphs. 

Unfortunately, Eqs. ( |2.4] ) do not possess any finite energy solutions. This is established 
by a variant of the argument relevant to the m? = case . 

Consider the symmetric tensor 

6'^ = 2 Tr (e'E^ + B'B^ - —{E^ + B^ + m^Al)] . (2.5) 



Using (|2.4| ) one verifies that for static fields dj 6^"^ = 0. Therefore 

/ d^r 6'' = I £r d^{x' e^') = I dS^x'e^' . (2.6) 

Moreover, the energy of a massive gauge field (with no mass for the spatial components) 
can be written as 



8=ld^r{-i:r[E^ + B- + ^.{D^E^) + Tr (mA^ + ^)^ • (2.7a) 



The second trace in the integrand enforces the constraint ( |2.4a| ). Consequently, on the 



constrained surface the energy is a sum of positive terms |]T0[ : 

g = Jd\ {- Tr [e' + B^ + m^Al) } (2.7b) 

and E, B and Aq must decrease at large distances sufficiently rapidly so that each of them is 



square integrable. This in turn ensures that the surface integral at infinity in (|2.6|) vanishes, 
so that static solutions require 



d^rO'' = . (2.8a) 

On the other hand, from (p.5|) , we see that 6^^ is a sum of positive terms 

r = - Tr {E^ + B'^ + 3m^Al) , (2.8b) 



hence (p.8|) imply the vanishing of E, B and Aq. 

The absence of finite energy static solutions can also be understood from the differential 
equations ( |2.4| ). Eq. ( |2.4aj ) possesses solutions for Aq that are either exponentially increasing 
or decreasing at infinity. Rejecting the former removes the freedom of imposing further 
conditions at the origin, and necessarily the exponentially damped solution devolves into 
one that is singular (not integrable) at the origin; see the Appendix. (This situation can 



be contrasted with, e.g., the magnetic dyon solution |jTT]], where absence of the mass term 
allows solutions for Aq with unconstrained large-r behavior, leaving the freedom to select 
the solution that is regular at the origin.) 

A similar argument shows that there are no "static" instanton solutions. These would 
be solutions for which t is replaced by —ix^, Aq by iA^ and presumably one would seek 
solutions periodic in 0:4 with period (3 = -^ = —J '^ J"' . An X4-independent solution is 



necessarily periodic; it would satisfy ( |2.4| ) with A/^ replacing Aq and opposite sign in the 
right side of (|2.4b| ). But analysis similar to the above shows that finite-action solutions do 
not exist. 



III. HARD THERMAL LOOPS FROM THE COMPOSITE EFFECTIVE ACTION 

In this Section, we present a derivation of the non-Abehan Kubo equation ( |1.1D based 
on the composite effective action of 0, a generahzation of the usual effective action (ob- 
tained by couphng local sources to the fields) in which one additionally introduces bilocal 
sources. In the QCD case, the composite effective action is given by S{A)+Tc{A, G^), where 
G(j){x,y) are (undetermined) two-point functions, and the labels (p = A,i/jX denote either 
gluons, or fermions-antifermions, or ghosts- ant ighosts, respectively (in the end, ghosts play 
no dynamical role, beyond maintaining gauge covariance of the final result). S{A) is the 
pure Yang-Mills action, and 

r,(A,G^) = ^( Tr InG^i + Tr V^'Ga^ 

-i( Tr InGl^ + Tr V:^^G^ + Tr lnG7^ + Tr Pr'^c) (3-1) 



when 2PI contributions are omitted. The trace is over space-time arguments as well as over 
Lorentz and group indices. The gauge coupling constant g, which was previously scaled to 
unity, is here reinserted. V^^ is computed from the usual QCD action Sqcd (e.g. in the 
Feynman gauge): 

The fields carry group and space-time indices, which are symbolically subsumed into the 
space-time labels x,y. 

The truncated composite effective action (|3.1| ) comprises the first, dynamical approxima- 
tion that we make and reflects the known fact that hard thermal loops arise from one-loop 
graphs. The full composite effective action of course coincides with the ordinary effective 
action when the two-point functions are evaluated by imposing stationarity requirements, 
and the above truncation reproduces the standard one-loop action involving Tr InPT^. 
Nevertheless subsequent analysis is more transparently organized in the composite effective 
action formalism. 



As indicated in [|^ , 5* + Fc is stationary for physical processes. This yields the conditions 

D^F^f" = J" , (3.3a) 

G-^^ = V-^\ = A^,C- (3.3b) 

Computing the local induced current J^{x) = —j^^ involves differentiating T)7^ with 
respect to A^. Since the P7^ depend locally on A, the resulting current is the local limit of 
a bilocal expression constructed from the two-point functions G^{x,y): 

J^(x) = limJ'^(x,i/) , (3.4) 

where the bilocal current J'^{x,y) = T"-Jl^{x,y) is given by 

nx,y) = g {T^^^^D:G/^x,y) + d^G^{x,y)) + zgr^tT^T'^G^x^y) (3.5) 

with r^^^ = 2g'^ga.y — gaQp-y ~ Q^Qpa- The trace "tr" is taken over Dirac spinor as well 
as internal symmetry indices, and we have defined GA,c,{^,y) = [r°'iT^]GA,Cab{,x,y) with 
D,GA{x,y) = d^GA{x,y) + g[[A{x),T%T']GAbc{x,y). 

We now use eqs. ( ^.31 ) - ( p3.5D to study "soft" plasma excitations. "Soft" means that both 
the energy and the momentum carried by a particle are of order gT, for a coupling constant 
(7^1, while particles with energy or momentum of order T are called hard (see e.g. 0). 
The strategy is to solve the system of coupled equations (|3.3| ), in order to derive from ( p.5|) 



the expression ( |1.3| ) for the local current J^. We approximate eqs. (|3.3| ) by expanding them 
in powers of g. The approximation scheme we use was first proposed in P] for deriving hard 
thermal loops from the Schwinger-Dyson equations. It represents an essential step in that 
derivation. Earlier work on the QCD plasma (in which this approximation was not used) 



is reviewed in |]T2[|. Following p[, we introduce relative and center of mass coordinates, 
s = X — y and X = ^{x + y), respectively. In these new variables the partial derivatives 
carry different dependences on g: dg ^ T and dx ~ gT. This comes from the fact that ds 
corresponds to hard loop momenta, whereas dx is related to soft external momenta. See |^ 
for a detailed account. 



Next, motivated by the expression ( |3.4| )- (|3l5| ) for the current, we expand G^f, in powers 



oig: 

G^ = Gf+gGf+g'Gf + ..., (3.6) 

where G^ is just the free propagator at temperature T and G^\ i > 1 are determined 
by (|3.3b|) . At leading order in g (to which we restrict ourselves in the sequel), the bilocal 
current (|3.5| ) depends on G^ and G^ : 

d^G^T^X, s) + r'A'^,{X)GT:^{X, s)) - d^GfliX, s) 



J!i{X,s)=g'f 



2 rabc 



■■- 0/37 



+tg'tTYT''G^2\x,s) + 6J!:{X,s) 



(3.7) 



where G(^(X, s) = G^{X + |,X — |) [and similarly for J{X, s)]. We have added the term 
6Jl^{X, s) in order to compensate for the loss of gauge covariance due to non- locality: 



5J!iiX,s)=g's-A\X) 



race rb 



pace fbcd I o f\s/^{0) A*!^, 



s) + d^Gf''\s)) + ^trr^TVG'J'^(s) 



. (3.8) 



Note that this term vanishes in the local limit. 

Now, we derive from ( [j.3b| ) a condition on Go, . [It turns out to be convenient to expand, 
instead of ( |3.3b| ), the equivalent equations V'^^Gfj, = G^j,!)^^ = /, in which we disregard 
temperature-independent contributions.] The 0((y')-condition does not fix G^ uniquely; 
hence we need to go to 0{g'^). The condition so obtained on G^ can be used to derive a 



constraint on the bilocal current. The subsequent derivation of this constraint [eq. ( p.l3[ )] 
is similar to the one given in |@|, to which we refer the reader for details. Momentum space 
is most convenient, i.e. 



G^{X,k) = I dSe'^-'G^{X,s) 



(3.9) 



the explicit forms for the thermal parts of the free propagators being {e.g. in Feynman 
gauge): 



G 



Aab 



Gf^'^'ik) 



Q{0)ab 



(A;) = -27r5V5(A;2)nB(A;o) 
{k)= 27T6''^6{k^)nB{ko) , 



(3.10a) 
(3.10b) 
(3.10c) 
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where ns^Fiko) = l/(e'^l'^°l =F 1) are the bosonic and fermionic probabihty distributions. 
Similarly, for the bilocal current in momentum space one writes 

J^'iX,k)= Jd^se"'-'J>'{X,s) . (3.11) 

In the limit s — ;> 0, or equivalently y ^ x, where X = x, 

nx) = r{X) = / ^ nX, k) . (3.12) 

The resulting constraint on the bilocal current is |]^: 

Q ■ Dx J^{X,k) = A7ig''Q^^QPk^Fp^6{k^)^[N nB{k^) + NpUFik^)] , (3.13) 

where Q'^ = g = (l,Q). 

Our next task is to make contact between ( p.l3| ) and the gauge invariance condition 
for the generating functional of hard thermal loops. Our strategy consists in transforming 
( p.l3| ) into two distinct conditions for positive and negative fco's, and in taking advantage 
of the symmetry properties that arise. We first integrate the equation ( |3.13| ) over |k| and 
fco > 0. Due to the 5{k'^) on the right side, the bilocal current is non-vanishing only when 
k^ = |k|; hence Q can be replaced by a unit vector g = t^. The integration thus yields: 

g+ • Dx JUX, q) = -2v^7r3m2g^g^F,o , (3.14) 

where we have defined 

//•CxD 
|kprf|k| / dkoJ>'{X,k). (3.15) 

Similarly, upon introducing 

J'^{X,q)= f\k\^d\k\ f dkoJ^'{X,k), (3.16) 

J J — CO 

the integration of ( p.l3| ) over |k| and k^ < gives: 

g_ ■ Dx J^iX, q) = -2V2Tx^m''Q^_QP_Fpo , (3.17) 

wherefrom one sees that J^{X, — g) satisfies the same equation (|3.14|) as J7+(X, g). 



Now, using J d^k = J dQ\\i\'^d\h\dko, we rewrite the expression (|3.12|) for the local current 
as J^{X) = f 72^[<^+(-^; q) + J^{X., q)]. Here, q can be replaced by — g in each term of the 
integrand separately, since q spans the whole solid angle. Therefore, we can write 

dq 



J^(X) 



(2^)^ 



J'{X,q) 



where J'^{X,q) is defined as 



(3.18) 



J^iX,q)^JUX,q) + J'^iX,-q) 



(3.19) 



and satisfies, as a consequence of (|3.14|) and ( |3.17|) 



g+ ■ Dx J^{X, q) = -4:V2r^m^Q^+Q'+F^ 



pO 



(3.20) 



From this, after decomposing 



J^(X, q) = J^{X, q) - AV^TX^m^Q'iA^ 



(3.21) 



we get as our final condition on the bilocal current: 

Q+ ■ Dx J^{X, q) = 4:V2r^m^Qld^xiQ+ " ^) • 

Let us now assume that J'^{X, q) can be obtained from a functional Vr(A, q) as 

5W{AA) 



(3.22) 



J'{XA) 



5A^{X) 



(3.23) 



Equation (|3.22|) then implies that W{A,q) depends only on A_|_, i.e. W{A,q) = W{A 
and J'^ = —sA^Q+- ^^ turn, 14^(A+) satisfies, as a consequence of ( p.22|) , 

6W{A+ 



Q+-Dx 



6A, 



4:V2 7r^m^dlA~ . 



x^ 



(3.24) 



By introducing new coordinates (a;+, a;_,x_L), 



x+ = Q- ■ X, z_ = Q+ ■ X, xx ■ g = , 



(3.25) 



we can rewrite Q+ ■ dx as 9+ and (|3.24|) becomes 
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\ 51V(^ 



bA, 



^^PlT^^m^d\A^ . (3.26) 



This equation was first derived in |0| , as an expression of gauge invariance of tfie generating 
functional for hard thermal loops, and has since then been studied by several authors. Here, 
it is seen to be a consequence of the stationarity requirement on the composite effective 
action. 

It has been shown in ||l| that iy(A+) is given by the eikonal of a Chern-Simons gauge 
theory. This observation is our last step towards deriving the approximate expression for 
the local current J^{x) in eq. (|3.3a|) . The subsequent development follows p and the result 



is exactly the non-Abelian Kubo equation (|1.1| ) with the form (|1.3| ) for the induced current. 



IV. CONCLUSIONS 



The behavior of the quark-gluon plasma at high temperature is described by the non- 



Abelian Kubo equation (IJ.) - (L3). We have studied the static response of such a plasma 
and proved that there are no hard thermal solitons (this result is supported and illustrated 
by numerical integration). The absence of "static" instantons is established by invoking 
a similar argument. In addition, the static non-Abelian Kubo equation indicates that the 



non-Abelian electric field is screened by a gauge invariant Debye mass m = gT J ^3^ . 

Furthermore, we have derived the non-Abelian Kubo equation from the composite ef- 
fective action formalism. Indeed, the requirement that the composite effective action be 
stationary leads, within a kinematical approximation scheme taken at the leading order, to 
the equation obtained in p[ by imposing gauge invariance on the generating funtional of 
hard thermal loops. 

Let us mention some problems deserving further investigation. Finding non-static solu- 
tions to the non-Abelian Kubo equation is an appealing — if difficult — task, since such 
solutions would correspond to collective excitations of the quark-gluon plasma at high tem- 
perature. Also, it would be interesting to investigate the next-to-leading order effects in 
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the kineinatical approximation and to see if they are gauge invariant; we hope that our 
formahsm is well suited for such an investigation. Furthermore, it is clear that rc{A,G^), 
when evaluated on the solution for G(j, obtained from ( |3.3b| ) and (|3.6|), coincides with the 
T{A) constructed from the Chern-Simons eikonal. While our derivation establishes this fact 
indirectly, an explicit evaluation of the relevant functional determinants in the hard thermal 
limit would be welcome. 

NOTE ADDED 



We have now seen recent papers [|T^ wherein the response equations are also analyzed. 
Moreover, local equations are found for time-dependent, but space- independent gauge fields, 
and for non-Abelian plane waves. The starting point of these investigations is a non-local 
expression for the induced current (see PJT3|]), 

-^v,j^ duUab{x,x-vu)wY!'{x-vu) , (Nl) 

which appears different from our local, but coupled, form ( |1.3|) - (|1.6| ). Here we exhibit the 
steps that explicitly relate the two. 

Beginning with our form for the induced current, ( |1.3|) - (|1.6| ), we observe that, owing 
to the integration over the angles of g, we may collapse these expressions into 

'^f{x) = m^ j^Q%[a^{x)-A^{x)) , (N2) 

where 

9+a_ + [A+, a_] = d-A+ . (N3) 

Eq. ( [IN 3D may be integrated, yielding 

("OO 

a'^(x) = / duUab{x,x-Q+u)d-A\{x-Q+u) . (N4) 

Here Uab satisfies 

12 



d 

— Uab{x,X -Q+U) = Uac{x,X - Q+u) fcbd ^'^{x - Q+u) 
Uabix,x) = 6ab ■ 



(N5) 



Also A°i may be presented as 



f 

Jo 



(iw -— <j Uab(x,x — Q+u) A''_(x — Q+u) 
du ' 



du Uabix, X — Q+u) < d+A^_ {x — Q+u) 
- f'^Alix - Q+u) A'lix - Q+u)}. 



(N6) 



[We have assumed that no contributions arise at infinity.] From ( |N2| ), ( |N4| ) and ( [N 6| ) , it 
follows that the induced current can be written as 



771 



dq 



— - j^(x) = mM -— QM du Uab{x, x - Q+u) Ft+{x - Q+u) 
I J 47r JO 



(N7) 



which coincides with the expression ( [Nil) derived in |T^, after the notational replacements 
m — > v^i^p, dq -^ dQ, (5+ -^ v^ and F_+ ^ v ■ E are performed. 

The time- dependent, space-independent equation found in [|13i is easily derived in our 



formalism, also. When there is no space dependence, eqs. (|1.6|) can be written as 



d+{a^ - A+) + [A±, a^ - A±] = 



(N8) 



and are solved by a^ = A+. Hence: 



'^r = ^l§iQ^-Q-nQ^-Q-rA. 



(N9) 



of which only the spatial component is non- vanishing: 



m 



■J / 4^ y y J 3 



(NIO) 



This coincides with the result in 13 



Similarly, the induced current for the non-Abelian plane wave in |T3[ corresponds to: 

Q± -p 



a+ 



Q+ -p 



A+ip-x) 



(Nil) 
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in our formalism, with p = {uj, k) being the corresponding wave vector. 

APPENDIX 

In this Appendix we analyze in greater detail and integrate numerically the radially sym- 
metric version of the static response equations (p.4|) , in the SU{2) case. Radially symmetric 
SU{2) gauge potentials take the forms: 



A"=r-r-r-)M^ + .-^F ^~'^^^^^ 






Al = f''^-^ , (Al) 

where a residual gauge freedom has been used to eliminate a term proportional to f'^f'^. 
We substitute the Ansatz (|A1| ) into ( p.4| ). The resulting equations give us the freedom 



to set one of the two (/),'s to zero; we obtain, 



x'^K = {K'-j'-l)K, (A2) 

where we have set 02 to zero, rescaled x = mr and defined J{x) = g{r), K{x) = 0i(r). 

We now investigate this system of coupled second-order differential equations. First, we 
see that they possess the following two exact solutions: 

J = 0, K = ±1 , (A3a) 

J = Joe"^ K = . (A3b) 



Eq. ( A3a ) corresponds to the Yang-Mills vacuum, while ( |A3b ) is the celebrated Wu-Yang 



monopole plus a screened electric field. 

In the asymptotic region a; — ;► cxd, the regular solution of the system (^) tends to (|A3aD 



with J approaching its asymptote exponentially. (Of course there is also the solution with 
J growing exponentially, which we do not consider.) 
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Near the origin, J and K behave either hke the vacuum ([A3 a]) or approach the monopole 
solution ( [A3b| ) as follows, 



J{x) —^Jq + ... , 

'277, X 

\ — 



/ZTT X \ 

K(x) -^ Koy/xcos[ — In— ]+..., (A4) 

V r XqJ 

where r is correlated with Jq as 



r = . . (A5) 

\/4 J2 + 3 

Only the vacuum alternative at the origin leads to finite energy. However, since we must 
choose one of two possible solutions at infinity (obviously we pick the regular one), the 
behavior at the origin is determined and can be exhibited explicitly by integrating the 
equations ( |A2|) numerically. Starting with regular boundary conditions at infinity, we find 
the profiles presented in Figure 1. They show that the monopole solution ( [A3b| ) is reached 
at the origin, with K vanishing as in ( [A4| ) - ( |A5| ), a result consistent with our analytic proof 
that there are no finite energy static solutions in hard thermal gauge theories. 
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